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EXAMINATION OF BOUNDARY CONDITIONS FOR HEAT 
TRANSFER THROUGH A POROUS WALL 
by Raymond S. Colladayand Francis S. Stepka 
Lewis Research Center 

SUMMARY 

An analytical model describing the heat transfer through a porous wall is presented. 
It utilizes an experimental upstream wall heat-transfer coefficient as a boundary condi- 
tion. An experiment was conducted to determine this coefficient for a surface simulating 
the surface of a typical full-coverage -film-cooled turbine blade, and the results are pre- 
sented. 

The temperature rises of the coolant approaching normal to the upstream surface 
and through a porous wall are more reasonable with the proposed model than with pre- 
vious models. 


INTRODUCTION 

An analytical model is developed in this report that can be used to predict the cool- 
ing air temperature rise and the metal temperatures through a porous wall. Previous 
analytical models for this purpose are presented in references 1 to 5. The most recent 
and comprehensive of the previous models (given in ref. 1), was utilized for the design 
of a full -coverage -film -cooled turbine blade at NASA. The results indicated that the 
rise in temperature of the cooling air as it approached normal to a porous turbine blade 
wall was greater than the rise through the wall even though the internal wall configura- 
tion was very effective from a heat -transfer standpoint. This apparent paradox led to 
the formulation of an alternative model which utilizes previously unavailable empirical 
data for the heat -transfer coefficient on the upstream surface of a porous wall. 

This report (1) presents the proposed model, (2) presents the results of a prelimin- 
ary experiment (conducted by the University of Minnesota under NASA contract) to deter- 
mine the heat -transfer coefficients for flow normally approaching a typical perforated 



wall configuration over a range of cooling air Reynolds numbers, and (3) compares the 
resulting model with that of reference 1. 

ANALYSIS 

In this section, a model for describing heat transfer from a porous wall is presented 
which utilizes, as a boundary condition, a Nusselt number correlation on the upstream 
face. Experimental data for this boundary condition were obtained by the University of 
Minnesota and are presented herein. 

The model in reference 1 is singled out for comparison with the presently proposed 
model because it represents the most recent and comprehensive treatment of the subject 
and did make an attempt, with some success, to overcome some of the shortcomings of 
earlier investigations (refs. 2 to 5) . 


Heat -Transfer Models 


Let model I refer to that presently proposed and model II to that of reference 1 . The 
energy balance and resulting equations for coolant and porous matrix temperature for 
region I (see fig. 1) are the same for both models. 

The heat-transfer process is described by assuming a one -dimensional model with 
constant properties where heat flowing by conduction through the wall in the negative 
x-direction is continuously transferred to the coolant in counterflow by convection. With 
the porous wall being composed of many small cooling passages within the matrix mater- 
ial, it is further assumed that conduction in the fluid in region I is negligible compared to 
that of the matrix. 

Consider an element in a porous wall of area A normal to the x-direction as shown in 
figure 1. An energy balance on this element requires that the net heat conducted by the 
matrix into the element must be equal to the heat transferred by convection to the fluid; 
hence 


d 2 T 


k e A 


w 


dx 


Ax = h m aA m< T w ' T c> 


where h m is the internal matrix surface -to -coolant heat transfer coefficient, AA m is 
the internal matrix surface area within the element, and k g is the effective thermal con- 
ductivity of the porous wall. (All symbols are defined in appendix A. ) Defining Z as 
the internal surface area per unit volume, 


Z = 


AA 


m 


AA X 


2 



leads to the following expression for the coolant temperature in region I: 


_ _ T k e d T w 

c " w ' k m z ** 

For a homogeneous pore structure, Z is constant with x. 

An energy balance on the coolant within the element in figure 1 gives 


( 1 ) 
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or 


k e dT w dT c 


G C ,2 dx 
c p dx 


( 2 ) 


Differentiating equation (1) with respect to x, neglecting entrance effects by assuming 
h m constant, and substituting into equation (2) leads to the following equation for the tem- 
perature distribution in the matrix material: 


d 3 T h Z d 2 T h Z dT 
w m w _ _m w 


= 0 


. 3 G C ,2 
dx c p dx 


k e dx 
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In dimensionless form, equations (1) and (3) become 


e = e - - 6" 

c w w 


(4) 
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where the prime denotes differentiation with respect to 4, and 
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(6c) 


(6d) 


(6e) 


where H m is the internal volumetric heat-transfer coefficient equal to h m Z . 

The solution of equations (4) and (5) leads to the following expressions for the coolant 
and wall temperature profiles in region I: 


where 



Boundary Conditions 


( 7 ) 

( 8 ) 

(9a) 

(9b) 


Model I. - Three boundary conditions are needed to evaluate the constants Cj, C 2 , 
and Cg. These constants are in turn functions of X and j3, which are known param- 
eters if the internal volumetric heat -transfer coefficient and effective thermal conduc- 
tivity are known for the porous wall configuration. 

The boundary conditions proposed are 
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( 11 ) 


e (l) = 1 

w v ' 

N0 W (O) = e'JO) (12) 

fl c (°) = ^ 0 w (O) (13) 

X 


where N = hX/kg. 

The first condition (eq. (11)) is that of a specified wall temperature T w Q . The re- 
maining two boundary conditions satisfy an energy balance at the interface at x = 0. It 
is assumed that the heat flux from the surface can be expressed as a product of a heat- 
transfer coefficient h^ and the temperature difference between the wall surface T w . 
and the coolant supply T . Continuity of heat flux then gives 


dT 


h i< T w,i- T e J 


w 


dx 


lx=0 


the dimensional form of equation (12). The heat-transfer coefficient must be obtained 
from an empirical correlation. This correlation should include the effect of coolant flow 
rate, geometrical hole array, hole spacing, hole diameter, etc. 

Even though the flow in the vicinity of a hole is three-dimensional, on a large scale 
(compared to the hole diameter and spacing) where average temperatures in the x = 0 
plane are desired, a one -dimensional analysis can be assumed. 

The heat loss from the x = 0 face can also be expressed as 


G C (T . - T ) = k 

C p V C,1 C,°°' € 


dT 


w 


dx 


|x=0 


which, when nondimensionalized, leads to equation (13). Note that the wall temperature 

T . and coolant temperature T . are a priori unknown. 
w> *■ ^ i 

The resulting expressions for the three constants in equations (7) and (8) are 

Cj = 0 (14a) 

N - a 9 

C 0 = — (14b) 

" ^ a 

(N - a 2 )e 1 - (N - a^e 2 
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(14c) 


N - a 1 

c - - 

'“'3 - 

cl< Qa 

(N - a 2 )e 1 - (N - a^e ^ 

It is customary to evaluate the thermal performance of a porous wall in terms of its 
effectiveness rj, where 


V = 


'c, o 


' w, o 


^c 

T c 


(15) 


Effectiveness is a measure of the efficiency of the porous wall as a mechanism for trans- 
ferring heat to the coolant. 

If the effectiveness is specified, then 

0 C U) = v (16) 

and the volumetric heat-transfer coefficient H m can be calculated as a function of flow 
rate G c by satisfying equation (16). Such a procedure, namely, that of satisfying a 
fourth boundary condition, is usually adopted in calculating the internal heat-transfer 
coefficient h . The internal matrix area per unit volume Z can often be found inde- 
pendently. 

Regarding the coolant temperature, note that only the total temperature rise from 
x = -o° to x = L is specified. How this total rise is distributed between regions I and II 
is a function of the model used to describe the heat -transfer process. 

With the constants of integration thus evaluated, equations (7) and (8) may be used to 
evaluate all fluid and matrix temperatures in the wall. Only the outer wall temperature 
(or heat flux to the wall), coolant supply temperature, effectiveness (or volumetric heat 
transfer coefficient), and effective thermal conductivity are specified. If a heat flux 
boundary condition is given at x = L with the wall temperature initially unknown, the 
boundary conditions for the model remain the same. However, to determine the dimen- 
sional temperature distributions, the exterior wall temperature must be calculated from 
the following overall energy balance equation: 

«"(!.) = G c C p ,(T w 0 -T c J 

If dimensionless temperatures are desired, only rj (or H ) and k g are required. 

Model H . - The boundary conditions for model II are presented in reference 1 and 
repeated in appendix B. Some of the symbols have been changed to be compatible with 
those of model I. 
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It should be noted that, in model II, the exponential coolant temperature profile for 
x < 0, given by equation (B3) occurs only if the flow in region II remains one -dimensional 
near the x = 0 face with no surface convection. However, an interface heat-transfer 
coefficient hu for model II can be defined as follows: 


« i ’=0 = h i 


T _ - T (0) 

c , CO W v ' 


= k 


^w, 

dx 


•x=0 


(17) 


or 


f e w(°) “ h i e w<°) 


Hence, 




+ C 3 a 2 

+ c 3 


(18) 


EXPERIMENTAL INVESTIGATION 
Apparatus 

The test section used to determine the heat-transfer coefficients for flow normally 
approaching a scaled -up model of a typical perforated wall of a turbine blade is shown 
schematically in figure 2. The surface of the wall at which the coefficient was deter- 
mined consisted of a 0.025-mm- (0.001-in.-) thick sheet of stainless -steel shim stock. 
This sheet was cemeted to a base plate of 1. 59-mm- (1/16-in. -) thick reinforced plastic, 
which in turn was bonded to a 7.62-cm- (3-in. -) thick layer of polyurethane foam insula- 
tion. This composite wall was a square 20. 3 cm (8 in.) on a side, which had 36 holes, 
each 6.35 mm (1/4 in.) in diameter drilled through the wall. The holes were arranged 
in a square array with 5-diameter spacing. Cemeted to the edges of the wall were 
1. 27-cm- (1/2-in. -) thick plywood sheets to form the sides of a plenum chamber for the 
entering air. 

The average heat-transfer coefficient at the surface of the perforated wall was ob- 
tained from the measured heat flux, the average surface temperatures, and the incoming 
air temperature. 
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The heat flux from the stainless -steel sheet (heater strip) was generated by electric 
current. This current was generated by a 6-volt battery. Leads from the battery were 
attached to copper bus bars that were soldered to opposite edges of the heater strip. The 
current flow was regulated by a variable resistor and determined by measuring the volt- 
age drop across a shunt. The heat input to the heater strip was then determined from 
measured current flow and the resistance of the heater strip. The average temperature 
of the heater strip was obtained by arithmetically averaging the readings obtained from 
five iron-constantan thermocouples (shown in fig. 2) that were cememted to the back side 
of the heater strip. The thermocouples were electrically insulated from the heater strip 
by a thin layer of copper oxide. The inlet air temperature was measured by an iron- 
constantan thermocouple located in the center and about 20 hole diameters upstream of 
the test wall. 

Ambient air was drawn to the test wall by means of a centrifugal blower downstream 
of the test section. The velocity approaching the wall was obtained from the measured 
flow, the wall area, and the inlet air density. 


Procedure 

The first set of data was obtained with the temperature difference between the heater 
strip and approaching coolant kept approximately constant while the air Reynolds num- 
ber was changed. This was done to maintain constant fluid bouyant forces. In order to 
determine the sensitivity to these forces, another set of data was obtained wherein this 
temperature difference was varied at a constant Reynolds number. 


RESULTS AND DISCUSSION 

Experimentally obtained heat-transfer data for flow normally approaching a perfor- 
ated wall are first discussed, and these data are then used in the discussion of the ana- 
lytical models for heat transfer through porous walls. 


Experimental Data 

Nusselt number Nu for flow normally approaching a scaled-up model of a perfor- 

S 

ated wall surface of a typical full-coverage-film -cooled turbine blade is shown, in fig- 
ure 3, to be linear on a log -log scale with flow Reynolds number Re . The characteris- 

tic dimension used in Nu and Re was the spacing between holes . The figure shows 

s s 
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that with no flow Nu due to natural convection was 5. 4 and that with flow Nu in- 

s s 

creased to between 25 and 35 for Re of 500 and 2000, respectively. All the data used 

s 

in this figure were obtained with a constant temperature difference of about 11 K between 
the wall and entering air. 

The results of the test showing the effect of bouyant forces are shown in figure 4. 
For this test Re was maintained constant at 1650 and the temperature difference be- 

O 

tween the wall and the entering air was varied between 6 and 29 K. The results in the 
figure indicate that the effect on Nu was not large. 

Only a limited number of variables were included in the empirical correlation for 
the heat-transfer coefficient presented herein. Additional experimental investigations 
will be required to include the effect of variables such as geometrical hole array, hole 
spacing, and hole diameter. 


Analytical Models 

For the purpose of comparing results from the two models, a typical full -cover age - 
film-cooled turbine blade wall composed of many small, closely spaced surface holes 
connected by a network of intricate torturous internal passages was chosen. The follow- 


ing conditions were assumed in the analysis: 

Diameter of surface hole, d, mm 0.5 

Hole spacing, S, diam 5 

Effective thermal conductivity, k g , J/(sec)(m)(K) 5.2 

Wall thickness, L, mm 0.5 

Effectiveness, r/ 0.7 


Frequently, the effectiveness of a porous wall is assumed to be unity. However, for 
turbine blade applications where wall thickness is small and the amount of area in contact 
with the coolant is limited r\ is usually less than 1 (refs. 1 and 2). In reference 6, values 
of t) from 0. 6 to 0. 8 were assumed as typical of many blade wall configurations. Over a 
range of coolant flow rates of practical interest, it has been observed that typical full- 
coverage -film -cooled walls exhibit effectiveness values which do not vary significantly 
with flow rate. For this study a constant value of 0. 7 was assumed. 

Internal heat -transfer coefficient . - There exists a unique value of internal volumet- 
ric heat-transfer coefficient H m which satisfies equation (16) for a given mass velocity 
at a specified effectiveness. For each model, the H m distribution as a function of G c , 
which satisfies an effectiveness value of 0.7, is determined and presented in figure 5. 

The magnitude of the internal coefficient is significantly greater for model I than for 
model II, and the slope is less. The ratio of the slopes for model I to II is 0.74. At a 
G c value of 0. 35 g/(sec)(cm ), representative of that expected in many turbine applica- 
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tions, the H m ratio (and consequently the h m ratio) is 1.93. The specific value of the 
slope is a function of the variation of t] with G c . The slope is greater for a constant 
effectiveness than it would be if 77 decreased with increasing G . 

V 

Back -face heat - transfer coefficient. - The In distribution used in the boundary con- 
dition given by equation (13) for model I was obtained from the experimental results given 
in figure 3 with the appropriate hole diameter and spacing. 

Although a back-face heat -transfer coefficient is not required for model II, it is an 
implicit result of the model and can be calculated by equation (18). The resulting In 
distribution, given as a function of G c , is shown in figure 6 along with the experimental 
value used in model I. The coefficient calculated from model II is considerably greater 
than would be expected on the basis of the experimental results. For example, at a G £ 
value of 0.35 g/(sec)(cm^), the he ratio for model II to I is 4.6. 

Actually, model II implies no convection on the x = 0 surface. However, the coef- 
ficient h^ is defined by equation (17) to give the convective equivalent of the back -face 
heat flux that occurs only by conduction upstream within the coolant. 

o 

Temperature profi les. - For a mass velocity of 0.35 g/(sec)(cm ), the dimension- 
less temperature distributions of the matrix and coolant for each model are given in fig- 
ure 7. The difference in the coolant temperature distribution for the two models is con- 
sistent with the relative magnitudes of the heat-transfer coefficients h^ and H m be- 
tween the models. 

The treatment of the x = 0 boundary or interface conditions has a much greater 

effect on the coolant temperature than on the wall temperature . The assumption in 

model II using the asymptotic value for 0 M (O)/0' (0) (eq. (B9)) strongly affects the heat 

w w 

flux at x = 0. However, since the wall temperature gradiant at x = 0 must be less than 
at x = L, the effect of this assumption on the wall temperature is small. It should be 
noted that, though model II forces the asymptotic relation at x = 0 for all G c , model I 
reduces to the same case if the coolant flow rate is large. 

The coolant temperature distribution for model II predicts that 72 percent of the total 
temperature rise occurs upstream of the back face, as compared to 17. 2 percent for 
model I. It was this distribution between the coolant temperature rise approaching the 
wall and through the wall that caused the initial concern when model II was applied to a 
similar cooling situation. 

Often the internal heat-transfer performance of a wall configuration alone, exclusive 
of the upstream (x < 0) region, is desired. In that case, a useful parameter is the wall 
internal effectiveness defined as 


T - T . 

77 - c >° C > 1 

'w rr. 


■ w, o 


- T 
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When related to the overall wall effectiveness 77 , this equation becomes 

_ 1 - <y°> 

" w " 1 - « c «» 

In the present example, model II predicts an ? 7 W value of 0.4, while model I gives 0.66. 

CONCLUDING REMARKS 

The experimental heat-transfer performance on the upstream surface of a perfor- 
ated wall with normally approaching flow was obtained. This perforated surface simu- 
lated the surface of a typical full-coverage-film-cooled turbine blade. 

An analytical model describing the heat transfer through a porous wall was proposed. 
It utilizes experimental data for the upstream heat-transfer coefficient as a boundary 
condition. As a result, reasonable temperature rises of the coolant approaching normal 
to the upstream surface and through the porous wall were obtained. 

Because of the limited number of variables included in the empirical correlation for 
the upstream surface heat -transfer coefficient presented herein, it would be desirable 
if the experimental study could be extended to include the effect of variables, such as 
geometrical hole array, hole spacing, and hole diameter. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 7, 1971, 

720-03. 
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APPENDIX A 


SYMBOLS 


parameters defined by eqs . (9a) and (9b) 
specific heat 
constants of integration 
hole diameter 

surface averaged coolant mass velocity, pu 

porous matrix internal volumetric heat -transfer coefficient, hZ 
heat -transfer coefficient for x = 0 face 
heat-transfer coefficient for inside porous matrix 
thermal conductivity 
wall thickness 

dimensionless heat-transfer coefficient parameter, h^L/k e 
Nusselt number 
heat flux 

Reynolds number 

distance between holes 

temperature 

approach velocity 

spacial coordinate normal to wall 

porous matrix internal surface area per unit volume 

parameter defined by eq. (B2) 

parameter defined by eq. (6e) 


T - T 

overall effectiveness, c ’ 0 — - c r- °° - 

T - T 

w,o c,°° 


wall internal effectiveness, 


c, o 


c, 1 


T - T . 
w,o c,i 


C, °p 


T - T 
w,o c,°° 


dimensionless temperature, 



X 


parameter defined by eq. (6d) 
v kinematic viscosity 

£ dimensionless x variable, x/L 

p density 

Subscripts: 
c coolant 

e effective wall 

i location at x = 0 

o location at x = L 

s hole spacing 

w wall 

00 location at x = - 00 

II coolant in region II 
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APPENDIX B 


BOUNDARY CONDITIONS FOR MODEL II 


The model of reference 1 (designated model II) attempts to describe analytically the 
behavior of the coolant in region II (see fig. 1) which, with accompanying continuity con- 
ditions at the interface between the two regions, allows a complete closed -form analyti- 
cal solution to the problem without explicitly specifying any heat flux condition at the 
x = 0 face. 

An energy balance on a fluid element in region II yields the following convection- 
conduction counterflow differential equation: 


d 2 T 


II 


dT 


dx 


= G C 

2 c P dx 


II 


(Bl) 


which in dimensionless form becomes 


e ii - “ e h = 0 


(B2) 


where 


G c C p L 

a = — 

k 


e n- 


T nW- T c,~ 
T w, o' 're- 


integrating equation (B2) and repeating equations (7) and (8) give for model II 


® n (?) = c 4 + c 5 e 




a< £ a 9 | 

V0 " C ! + C 2 e 1 + C 3 e 2 


(B3) 

(?) 
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e c tt) = c 


1 



(8) 


Five boundary conditions are needed to solve for the constants of integration, and 
there is an added sixth condition if the volumetric heat-transfer coefficient is not speci- 
fied. The five boundary conditions are 


e (i) = l 

(ID 

e n (-«) = o 

(B5) 

e n (o) = e c (o) 

(B6) 

6' (0) = — 0{ T (O) 
w a/3 11 

(B7) 

0' c (0) = -A- 0*j( 0) 

(B8) 


a/3 


The first two boundary conditions are for specified outer wall temperature and cool- 
ant supply temperature, respectively. Equations (B6) and (B7) are x = 0 interface con- 
ditions requiring continuity of coolant temperature and heat flux, respectively. In equa- 
tion (B7), note that k c /k e = X/a/3. 

The interface condition in equation (B8) requires more explanation. Because con- 
duction through the coolant in region I was neglected, a discontinuity in the slope of the 
coolant temperature at x = 0 will result. A relation between these two slopes ( 0^(0 ) and 
0 jj('O)) can be found by combining equations (2), (Bl), and (B7): 


dT_ 


dx 


x=0 


G C 
c p 



(B9) 


It is then assumed that the ratio of second to first derivatives of the wall temperature at 
x = 0 can be obtained from the limiting case where the coolant and matrix temperatures 
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in region I are everywhere equal, that is, hr = 00 and r\ = 1. For this classical case 
(see ref. 3), equation (2) reduces to 

A; _ °c C p dT = 0 
dx 2 k e ^ 


where T = T = T, and the boundary conditions become 

W C 


dT 


dx 


= G C 
c p 

T(0) -T c oo 

x=0 



T(L) = T 


w 


The solution in dimensionless form is readily obtained as 




m = e 


- ^ 


(BIO) 


and 


0"(O) _ X _ ^w^ 
0’(o) /3 eyo) 


(Bll) 


Substituting equation (Bll) into (B9) yields the interface condition given in equation (B8), 
The resulting expressions for the five constants of integration are 

c t = 0 


C 2 “ 


b 2 (l - b 2 ) 


C 3 


cl cl 

b 2 (l - b 2 )e 1 - bj(l - bj)e 2 
b t (l - b t ) 


3. 2 

b 2 (l - b 2 )e 1 - bj(l - b 1 )e 


(B12) 
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where 


C 


4 


= 0 


C 5 ~ b l C 2 + b 2 C 3 


b 


1 


/3a x 
X 


b 


2 


X 


The constants Cg, Cg, and Cg are functions of /3 and X, which are in turn a 
function of H m> However, the local matrix and coolant temperatures given by equations 
(7), (8), and (B3) can be evaluated without explicitly specifying H m if the effectiveness 
condition (eq. (16)) is satisfied. 
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Figure 1. - Porous wall temperature profile model. 





Nusselt number, NU = h:S/K 



Figure 2. - Schematic of test section used to determine heat-transfer 
coefficients for flow normally approaching a perforated wall. 
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Figure 3. - Nusselt number for flow normally approaching perforated 
wall. Temperature difference between wall and coolant, ~11 K); at 
Reynolds number of 0, Nusselt number is 5.4. 
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